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O ■ Abstract 
O ; 

. A many-body Hamiltonian describing a system of Z protons and N neutrons moving in spherical 

^ : 

D , shell model mean field and interacting among themselves through proton-proton and neutron- 

' neutron pairing and a dipole-dipole proton-neutron interaction of both particle-hole and particle- 

psj . particle type, is treated within a fully renormalized (FR) pnQRPA approach. Two decoupling 

> : 

Tj" , schemes are formulated. One of them decouples the equations of motion of particle total number 

o . 

' conserving and non-conserving operators. One ends up with two very simple dispersion equations 

■ for phonon operators which are formally of Tamm-Dancoff types. For excitations in the (N-1,Z+1) 

system, Ikeda sum rule is fully satisfied provided the BCS equations are renormalized as well and 
\ therefore solved at a time with the FRpnQRPA equations. Next, one constructs two operators 

^ ■ ^1/^) '^i,-^t(— )^^^ which commutes with the particle total number conserving operators, and 

• • . )^~'^, and moreover could be renormalized so that they become bosons. Then, a phonon 

^ \ operator is built up as a linear combination of these four operators. The FRpnQRPA equations 

' are written down for this complex phonon operator and the ISR is calculated analytically. This 

formalism allows for an unified description of the dipole excitations in four neighboring nuclei (N- 
1,Z+1),(N+1,Z-1),(N-1,Z-1),(N-|-1,Z+1). The phonon vacuum describes the (N,Z) system ground 
state. 
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I. INTRODUCTION 



Understanding nuclear physics phenomena in terms of single particle motion in a mean 
field and nucleon-nucleon interaction is one of the most appealing aims of theoretical nuclear 
physics. Along the latest few decades many progresses have been obtained in the many body 
treatment of nuclear system. Among the chief achievements one distinguishes the Hartree- 
Fock, BCS and random phase approximation (RPA) 

HQ. 

The last mentioned approach 

can be formulated on the top of either a HF or a BCS ground state. The second version 
is conventionally called quasiparticle random phase approximation (QRPA). Among many 
successes of many body theories are the description of the ground state properties, of the 
excited states, including those of high spin of ground as well as of excited bands, of the 
electromagnetic transitions, and of the reaction mechanisms. 

The proton neutron interaction has been also treated. For example the proton-neutron 
pairing interaction is currently investigated for nuclei close to the drip line. It is still an 
open question whether there exists a T=0 proton-neutron pairing j^l phase, although there 
are some claims that the strong back-bending seen for ^^Fe could be explained by breaking 
a proton- neutron T=0 pair . 

The dipole proton-neutron interaction in the particle-hole (ph) as well as the particle- 
particle {pp) channels have been treated within the pnQRPA formalism in order to get 
a quantitative description of the double beta decay with two neutrinos in the final state. 
Of course, this phenomenon is very interesting by its own, but exhibits a special 
attraction for theoreticians having in view that it provides a test for the nuclear matrix 
elements which are also used in the neutrinoless double beta decay (Oz//5/3) calculations. In- 
deed, the discovery of Oz//3/5 process would answer a fundamental question whether neutrino 
is a Majorana or a Dirac particle. 

Standard pnQRPA calculations based on ph two body interaction yields a too large rate 
for the 2u(]i3 process. In Ref . |7] it was pointed out that the l3~^ transition matrix element is 
very sensitive to the pp interaction strength. Since the double beta decay has a branch which 
could be looked at as the hermitian conjugate matrix element of the virtual transition 
of the daughter nucleus to an dipole l"*" state in the intermediate odd-odd nucleus, many 
groups working on 2zy/9/3 decay included the pp two body dipole interaction in the pnQRPA 
calculations. By contrast to the ph interaction, which is repulsive, the pp interaction has an 
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attractive character. Since such an interaction is not considered in the mean field equations, 
the approach fails at a critical value of the interaction strength, gpp. Before this value is 
reached the Gamow- Teller transition amplitude, denoted by Mgt, is rapidly decreasing and 
after a short interval is becoming equal to zero. The experimental data for the amplitude 
Mgt is met for a value of Qpp close to that for which Mgt vanishes and moreover close 
to the critical value of gpp where pnQRPA breaks down. It is obvious that in this region 
of gpp, the pnQRPA results are not stable against adding interaction terms which are not 
encountered by the adopted many body approach. In order to restore the ground state 
stability, one needs to make a higher RPA calculation. The first formalism devoted to this 

n □ 

feature includes anharmonicities through the boson expansion technique [8|, |9[. Another 
method is the renormalized pnQRPA procedure (pnRQRPA) |Ji)j which keeps the harmonic 
picture but renormalizes the bosons by effects coming from some terms of the commutator 
algebra which are not taken into account in the standard pnQRPA. 



In a previous paper we have shown that the pnRQRPA does not include the higher 
RPA effects in a consistent way. Indeed, in Ref. [l^ only the two quasiparticle dipole opera- 
tors are renormalized due to the non- vanishing ave rag e of the quasiparticle number operators 
in the renormalized QRPA ground state. In Ref.|ll| we showed that having non- vanishing 
proton and neutron quasiparticle numbers, the scattering terms are also renormalized so 
that they finally satisfy bosonic commutation relations. Thus, a new phonon operator could 
be defined which includes, in addition, the scattering terms. We baptized the new renor- 
malization formalism as a fully renormalized proton-neutron quasiparticle random phase 
approximation (FRpnQRPA). The FRpnQRPA equations determining the phonon operator 
amplitudes have twice as much solutions, 2Ns, as the standard pnRQRPA equation. For 
half of these solutions one expects that the amplitude of the scattering terms prevails. In the 
quoted reference we pointed out which contribution is coming from the ph and which comes 
from the pp interactions in building up the new modes. It was shown that for charge conserv- 
ing excitations these modes are spuriousjl^. In a later publication we treated the scattering 
terms semi-classically and the harmonic mode defined there describes the wobbling motion 
of the isospin degrees of freedom jisl. 

It is worth mentioning that both renormalization approaches, pnRQRPA and FRpn- 
QRPA, violates the Ikeda sum rule (ISR)Q], saying that the difference of the total (3~ and 
(3~^ strengths for a single beta minus emitter is equal to 3(N-Z). Aiming at conciliating the 
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two features, renormalization and having ISR satisfied, recently a new phonon operator 
was defined whicli includes also the scattering terms but commutes with the particle total 
number operator. The resulting formalism is called fully renormalized pnQRPA,i.e. by the 
name that we adopted four years before. A distinct way to renormalize the RPA equations 
without modifying the number of solutions and moreover without deriving an explicit ex- 
pression for the phonon operator was formulated in Ref . jiojl . The method is known under 
the name SCRPA (self consistent RPA). 

n 

The present paper continues the investigation we started in Ref jll| , addressing the follow- 
ing issues. We want to see to what extent it is possible, in a quasiparticle BCS framework, 
to separate a harmonic boson operator conserving the particle total number from the com- 
ponents which do not commute with the total number operator. We shall see that the 
phonon operator conserving the particle total number, treated in a restricted A^^ dimen- 
sional space, can be embedded into an operator space acting on a 2Ns dimensional space 
and moreover has exactly the same properties as the phonon operator we have previously 
introduced in Ref.^^. Also, we aim at treating on an equal footing the particle-hole and 
deuteron like dipole excitations. For both cases the FRpnQRPA equations are written down 
explicitly. Analytical expressions for the renormalization factors will be derived. An impor- 
tant question addressed in this paper is whether having a particle total number projected 
phonon operator assures that the ISR is fulfilled. Another problem treated in the present 
paper is the following one. In the FRpnQRPA approach, the phonon operator is a lin- 
ear combination of the renormalized operators (see the notations from the next sections) 
A[^{pn),Ai^_^{pn){—Y^^,B\^{pn),Bi^^^{pn){—Y~^, neglecting the commutators of A and 
B operators. Here we define two renormalizable operators TZ\^{pn) ^TZi,^f^{pn){—Y^^ which 
are exactly commuting with the operators A\^{pn) and Ai^^^{pn)[—Y~^^ , where the latter 
operators are preserving the total number of particles. The phonon operator is defined 
as a linear combination of the A and TZ operators and explicit equations for the phonon 
amplitudes and energies will be derived. 

The objectives sketched above will be accomplished according to the following plan. In 
Section 2 we briefly review the main ingredients of the FRpnQRPA formulated in Ref. fll | . 
In Section 3 a consistent decoupling scheme is defined which results in providing two inde- 
pendent sets of Tamm-Dancoff equations describing excitations of the (N,Z) system in the 
nuclei (N-1,Z+1) and (N+1,Z+1), respectively. For the first type of excitation the Ikeda 



4 



sum rule is fulfilled. In Section 4, a complex phonon operator is defined, describing in a 
unified fashion states in the systems (N-1,Z+1), (N+1,Z-1), (N+1,Z+1), (N-1,Z-1), respec- 
tively. Finally, in Section 5 the important results are summarized and some conclusions are 
drawn. 



II. BRIEF REVIEW OF THE FULLY RENORMALIZED PNQRPA FORMALISM 

The single and double beta Gamow-Teller transitions can be described by the following 
many body model Hamiltonian: 



T;jm T;j,j' 

+ 2x E /?;W/3^m(A')(-)'^-2xi E prM)Pi%ip'^')i-r- (2.1) 

pn,p'n';ix pn,p'n';n 

where cl-j^ {cTjm) denotes the creation (annihilation) operator for a nucleon of type r (=p,n) 
in a spherical shell model \nljm). The time reversed state corresponding to \nljm) is 
\nljm) — \nlj — m)(— )^~"^. For the sake of simplifying the notation here only the quantum 
numbers j, m are specified. Also, the following notations have been used: 



'^pjmpnj'm' 

E(/>i^k,.|^jW)4„,ct^, = ^4-{pj\\a\\nj') 



1/i 



Vl -I 



(2.2) 



In what follows, sometimes the notations are simplified even more and the set of quantum 
numbers for a single r particle state will be denoted by r, the other quantum numbers being 
mentioned only if necessary. 

The Hamiltonian H is treated first by the BCS formalism which defines the quasiparticle 
representation through the Bogohubov-Valatin (BV) transformation: 

Ul + Vl = 1, sjm^i-y-"', r^p,n (2.3) 
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The images of the operators and Pf^ due to the BV transformations are: 



P:{k) 



(2.4) 



where the operators A\A,B\B , are the two quasiparticle and quasiparticle density dipole 
operators, respectively: 



The factors ^, ^, 7777 have the expressions: 



(2.5) 



Vk 



Jp 

i 
i 



{jp\W\\jn)Uj^Vj„, 
{jp\W\\jn)Uj^Mjn, 



Jp 

i 

Jp 
i 



Opl|o^llin)1^jpf/i„, J 
{jp\W\\jn)Vj^Vj„. 



2j + 1, 



(2.6) 



It is worth mentioning that throughout this paper the Rose convention for the reduced 
matrix element is adopted ^|. Within the proton-neutron quasiparticle random phase 
approximation (pnQRPA), one assumes that the two quasiparticle dipole operators ^41^, Ai^i 
satisfy quasi-bosonic commutation relations, while the quasiparticle density dipole operator 
commute with each others. Also the A operators commute with any of B'^ and B operators. 
The renormalized pnQRPA approximate the commutator Ai^,A\^ by keeping from its 
exact expression only the C number term and the monopole term. Moreover, replacing 
the adopted expression for the commutator by its average on the pnQRPA ground state, 
unknown yet, the operators A and A'^ can be renormalized so that the new operators obey 
bosonic commutation relation. Note that other commutator equations, are kept as in the 
standard pnQRPA, i.e. they are taken equal to zero. 

The idea advanced in Refs. |Q was that also the proton-neutron quasiparticle density 
dipole operators can be renormalized so that finally they are gifted with boson properties. 
Therefore, the mutual commutation relations assumed in Ref.jll[ 



are: 



Ai^{k),A\^,{k')] ^ 6k,k'Sl^,l^' 



1 - 



N„ _Np 

Jp 
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Blik),Al,{k') 
B,,{k),BUk') 



Jn 



Jp 



k — {jp, jn)- 



(2.7) 



Let us denote by Cl^\ and C)^'^ the averages of the r.h. s. of the first and third commuta- 
tion relations ()2.7|) on the pnQRPA vacuum state (|0)), respectively. Then, the renormalized 
operators are: 



.(2) 



Alik) 



Blik) 



=Al{k), A,,{k) = -j^A,,{k) 



(1) 



Blik), B,,{k) 



-Bi^{k). 



(2.8) 



Recalling the specific RPA convention to replace the operator commutators by their average 
on the ground state, one readily obtains that the renormalized operators defined above 
satisfy boson like commutation relations: 



'A^,{k),A\^,{k' 
B,,{k),BUk' 



Sk,k'Sti,tj.' fk, fk 



sign{C, 



(2)> 



(2.9) 



Further, the renormalized operators are used in order to define the pnQRPA phonon operator 



rl = Y: [xik)Alik) + Z{k)Dl{k) - Y{k)A,^,{k){~f-^ - W{k)D,,,{k){- 



(2.10) 

where D\^{k) stands for either B\^{k) or Bi^{k) depending on whether is equal to +1 or 
— 1. The amplitudes X, Z, W are determined by the fully FRpnQRPA equations provided 
by the operator equations: 



r 



r 
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(2.11) 



The number of FRpnQRPA equations is double the number of standard pnQRPA equations. 
They have to be solved at a time with the equations defining the constants C^^-* and C^^-*. 
The phonon vacuum is the ground state of the system which has the property that the 
corresponding average of the quasiparticle number operator is non-vanishing. Solving these 
equations, one finds out that half of the solutions have the amplitude X dominant while 
for the remaining ones, the amplitudes Z prevail. The latter solutions have been separately 
studied in Ref.[l3|. It has been proven that such solutions describe a wobbling motion of 
the isospin degrees of freedom. 



III. RESTORING THE GAUGE SYMMETRY FOR THE FULLY RENORMAL- 
IZED PNQRPA SOLUTIONS 



The major component of vacuum state |0) is a state characterizing the even-even system 
(N,Z) under consideration. Considering the inverse BV transformation for the dipole oper- 
ators A^,A,B\B one can easily check that one phonon states are mixtures of components 
describing the neighboring nuclei (N-1,Z+1),(N+1,Z-1),(N+1,Z+1),(N-1,Z-1). The first two 
components preserve the total number of nucleons but violate the third isospin component T3 
while that remaining components violate the total number of nucleons but preserve Ts.The 
latter two components mentioned above are the ones which are responsible for the violation 
of Ikeda sum rule (ISR), which is vahd for single beta decays of the (N,Z) nucleus. Aiming 
at getting a suitable structure for one phonon state, so that ISR is obeyed, it is desirable to 
start with linear combinations of the basic operators , A, B\ B which excite the (N,Z) nu- 
cleus to the nuclei (N-1,Z-M),(N-M,Z-1),(N-M,Z-M),(N-1,Z-1), respectively. One can check 
that such operators are: 

A[^{vn) = UpVnA{^{pn) + UnV^A^^.^{pn){- f-^ + UpUnB\^{pn) - V;,K5i,_^(pn)(-)i-^ 
Ai^iwi) = UpVnA.^ipn) + + UpUnBii^ipn) - VpVnBl_^{pn) 



Ai^(pn) = UpUnAi^{pn) -VpVnAl _^{pn){-] 



UpV^B,^{pn) - VpU^Bl_^{pn){-Y-^. 

(3.1) 



Indeed, expressed in terms of creation and annihilation particle operators the above operators 
are: 



, Aif,{pn) = - 



^p^n 



t 

1m' 



(3.2) 



In terms of the new dipole operators, the previously introduced one body operators become: 



P^tipri) = j{jp\Mjn)A\i,{pn), (3+{pn) = ■^{jp\\a\\jn)Ai,-^,{pn){-f, 



PrM) = f Oplkl|jn)Al^(pn), P;(pn) = |(j^||a||j„)Ai,_^(pn)(-)^ (3.3) 
The model Hamiltonian can be written as: 

H = ^rjalj^arjm + '^xYl (^pn;p'n'Al^{pn)Ai„{p'n') - 2xi J2 (Tpn■,p'n'A.l^l{pn)M^c{p'n'), 



JpJp' 



{jp\W\\jn){jp'\W\\jn')- 



(3.4) 



In order to study the harmonic modes which might be defined with the operators 

Ai^{pn) , Aifj_{pn) , A\^{pn) , Ai^{pn) we need their mutual commutation relations. These 
are given in Appendix A. Following the prescription of the fully renormalized pnQRPA, 
these commutators can be approximated as: 



Ai^{pn),AL{p'n') 



Sua' Si 7 , 



' " '■•2 P 



3n 

C/2 _ y2 _ 



Jp 

c/2 - y2 / 

a2 P 



Ai^,{pn),Al^,{p'n')] ^ 5M,M'^jp,v^jnj„,f^nK 
Alipn),Al_^,{p'n')i-y-^''] ^ 5^,^,6,^,, J,„,^,UpVp 



1-2- 



Jn 



1-2^ 
Jp 



(3.5) 



Here, iVp and Nn stand for the proton and neutron quasiparticle number operators, respec- 
tively: 

= X) alj^a^jm; r^p.n. (3.6) 



jm 



For what follows it is useful to introduce the notations: 



D,{pn) = - U'^ + ^j^iUl - K?)(iV„) - - ^P^^^^^^ 



Dp = UpVp 



1 - 



+ 1 

2 

2jp + r ^ 



(3.7) 



Here the average value of an operator O on the renormalized phonon operator vacuum state, 
undefined for the time being, is denoted by (O). According to Eq. (A.l) the quantities Dn 
and Dp represent to some extent the overlap of the ground state of the (N,Z) system with 
states describing the (N-|-2,Z) and (N,Z-|-2) systems, respectively. The latter states are 
obtained by adding one pair of neutrons and one pair of protons to the (Z,N) ground state, 
respectively. Both pairs, of neutrons and protons, have a vanishing total angular momentum. 
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Note that for the single particle states involved in single P" transitions, Di{pn) is a 
positive quantity, while D2{pn) might be either positive or negative, depending on the rel- 
ative values of the occupation probabilities for protons and neutrons in the states p and n, 
respectively. Let us denote by jc the critical angular momentum having the property: 



If Spj < epj^, then D2(pn) < 0, 
If epj > epj^, then D2(pn) > 0. 



(3.8) 



It is useful to introduce the functions: 



1 if D2(pn) > 0, 

if D2(pn) < 0, 

1 if D2(pn) < 0, 

if D2(pn) > 0. 
Then the operators ^1^, can be renormalized as: 



(3.9) 



\lDi{pn) 
1 



\D2{pn)\ 

One can certainly introduce new operators having proper boson like normalization: 



(3.10) 



Bi^{pn) = ep^A^^ + dp^A\A-f-^, 



(3.11) 



The result is that we have two pairs of operators satisfying boson like commutation relations: 



Ai^,{pn),A\^{p'n' 
B,,{pn),Bl{p'n' 



(3.12) 



In this Section we suppose that Dp and D„ arc negligible small. Also, we ignore the contribu- 
tion of the pairing interaction coupling the equations of motion associated to the operators 
A\ij_ and A|^. Under these circumstances the two boson operators are independent from 
each other, 

(3.13) 



Bi^(pn),A\^(p'n') 



0, 
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and consequently the equations of motion for the quoted operators are decoupled . 
'H,A,{pn)\ = [EpK - V^) + E^{V^ - UD] A{^{pn) + 2xY.4lv.nAM^i)^ (3-14) 

Here the following notations have been used: 



(7, 



[k) 

pn;p\n\ 



^(p||a||n)|D,(p,n)|V2^(pi||a||ni)|D,(pi,ni)r/', ^ = 1,2 



We look for the linear combination operators 



(3.15) 



(3.16) 



p,n 



with the properties: 



t ■ 

5 



t 



(3.17) 



These equations provide homogeneous systems of linear equations for the amplitudes X and 
Af, respectively. The compatibility condition yields for oj and VL the dispersion equations: 

^-2^[{p\\a\\n)fD,{pn) 



2xE 



2X1 E 



^J^^({p\\a\\n)fD2{vn) 



+ 1 = 0, 



1 



0. 



(3.18) 



The distinct feature of these equations consists of that the poles of the functions involved 
in the l.h.s. of Eq. ()3.18p are simple. Thus, while in the standard pnQRPA the dispersion 
equation depends on the squared energy, here the dependence is simply on energy. The 
important consequence is that the equations derived above do not have vanishing solutions. 

The states created by acting with the phonon operators and on the ground state of 
an (Z,N) system are dipole states l"*" in the neighboring odd-odd (Z+1,N-1) and (Z+1,N+1) 
nuclei. While the first state has a particle-hole (ph) character, the second one is a particle- 
particle (pp) like excitation. The phonon amplitudes have the expressions: 

^f{p\\(r\\n)Dy^{pn) 



X{pn) 
X{pn) 



-2x 



1 



+ i?„(V;2-[/2)-u; 

(T\\n)\D2ipn)\'^/^ 



Ep{U^ - V,') + E^{Ul - v^) - n 
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s. 



(3.19) 



The factors S and S are determined by the second equation from ()3.17p . requiring that the 
phonon operators are normahzed to unity. The final resuhs are: 



^ [{p\\a\\n)f D,{pn) 



1/2 



1/2 



(3.20) 



Note that both the dispersion equations and the phonon amphtudes depend on the renor- 
mahzation factors D, which at their turn depend on the quoted amphtudes. Therefore, the 
FRpnQRPA equations and the defining equation (3.7) should be solved at a time. In order 
to do that we have to provide a recipe for how to calculate the averages for the quasiparticle 
number operators involved in Eq. (3.7). In what follows, we shall describe separately the 
results for the phonon operators T"^ and As we proceeded in Ref.^|, we seek for a boson 
representation of the quasiparticle number operators, determined so that their commutation 
relations with the phonon operators are preserved. 



k,fj, 
k,fj, 

The expansion coefficients have the following expressions: 

Ck{r) = -(0| [[iV.,ri^(A;)l ,rI^(A;)l |0), T = p,n. 



(3.21) 



(3.22) 



Here |0) stands for the FRpnQRPA vacuum. Calculating the commutators involved in the 
above equations, one finds: 



k 

_ ^ \Xk{pn)\^ 
~ h Diipn) 
(iV„) = ECfc(n) 

k 

_ |Xfc(p?2)P 



^ p ' n ' ^ n ' p ' 



(Nn) 



a2 



p 



jp 



p n ' n p 



+ ^-^{ul-v^){ul-v;) 

Jn 



Jp 



.(3.23) 



These equations allow to express the averages {Np) and {Nn) as functions of Di{pn). Insert- 
ing the results in (3.7), a set of equations for the renormalization factors is obtained. Since 
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the averaged quasiparticle numbers depend on the X amphtudes, so do the renormahzation 
factors Di{pn). 

As for the deuteron hke phonon operator the results are: 

k,fi 

7V„ = J2Fk{n)gi^.{k)gl{k), (3.24) 

where 

Fk{r)^-{0\[[Nr,gi,.{k)],gl{k)]\0), r^p,n. (3.25) 
Using the final results for the coefficients F, one obtains: 

(TV,) = Y.Fk{p) 

k 

^ Do(vn) I'^p'^n^V'^n ;-2 '^'^n ^^n A^p ;-2 

(7V„) = X:Ffe(n) (3.26) 



k,p 



These can be viewed as a system of nonlinear equations for (iVp) and (iV„) for a given set of 
amplitudes X. Solving these equations and inserting the solutions in Eq.(3.7) one obtains 
the values of D2{pn). An easier way would be to express (iVp) and {Nn) as function of 
D2{pn) and then Eq.(3.7) becomes a nonlinear equation for D2{pn) that should be solved 
at a time with the equations (3.18) and (3.19). 

Consider, again, the equations associated to the operators and A. The quantities 
Di{pn) should fulfill an additional consistency equation caused by the fact that the phonon 
operator commutes with the quasiparticle total number operator. Requiring that also T3 
is preserved in the average, one obtains that the numbers of protons and neutrons are 
separately preserved. The average values of the proton and neutron number operators 
corresponding to the renormalized pnQRPA ground state have the expressions: 

[/2 _ y2 ^ 

p Jp 

N = Y.3l{V^ + ^^S^ W), (3.27) 

P Jn 
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These equations suggest that in order to have a pnQRPA ground state with non-vanishing 
quasiparticle numbers and with good proton and neutron numbers it is necessary to renor- 
mahze the BCS equations for chemical potentials. In this way the BCS and fully FRpnQRPA 

equations are coupled to each other. It should be mentioned that the above equations have 
been derived by assuming that the averages of the quasiparticle pair operators are vanishing. 
There is still an alternative option, namely to keep the standard BCS equations as they are 
but enforce the consistency restrictions: 

Ef^p^p3p((K4]oo + Mp]oo)) = E(^p 
p p 

E UnVnJn{{[aial]oo + Moo)) = E(Kz' ' U^Wn) : (3.28) 
n n 

As we shall see a bit later, following the second path sketched before, the ISR is violated. 

Despite the fact the phonon operator defined before is of a Tamm-Dancoff type, it might 
be viewed as a FRpnQRPA phonon operator acting in the extended 4Ns dimensional space. 
Therefore, the present results reproduce a set of A^^ solutions, with Ng denoting the number 
of proton-neutron states (p,n) which can couple at an angular momentum equal to unity, 
out of 2Ns positive solutions provided by FRpnQRPA. 

Indeed, denoting by 

the phonon operator can be expressed in terms of the renormalized operators: 

Alipn) = \ Alipn), Bl{pn) = ^ Slipn). (3.30) 

rL = E^(H^r'^'(H {u,VnDfA\^{pn) + Ury,DfA,,_,{pn){-f-^ (3.31) 
+ UMDB\^'^B\^{pn) - VMDB\^'^B^,-t.{pn){-f-^) 

= E [X(pn)4^(pn) - Y{pn)A,^^^{pn){-f-^ + Z{pn)B\^{pn) - ■W{pn)B,,_,{pn){-Y~^ 
The norm for the extended phonon is: 

E (X2(H - Y2(pri) + Z'^{pn) - W^{pn)) 
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p,n 



pn 



Di{pn) 
1 

Di{pn) 



p n n p 



p n 



n n 



(Nn) 


{Np) 


Jn 




(Nn) 


{K- 


Jn 





^ p^ n ' p ' n 





{Np)\ 







Jp 



Y,X\pn) = l. (3.32) 



p,n 



Thus, the Tamm-Dancoff phonon operator can be embedded in the operator space acting 
on the fully renormalized pnQRPA states. Reciprocally, Eq. (j3.32p tells us how to obtain the 
Tamm-Dancoff phonon amplitudes as functions of the fully renormalized pnQRPA phonon 
amplitudes. 

Note that, although the operator A'^ is a particle hole operator in particle representation, 
the Tamm-Dancoff equations accounts for quasiparticle correlations and therefore describe 
the system's small oscillations around a static BCS ground state. 

Let us check now whether the Ikeda sum rule is obeyed by the ground state corresponding 
to the Tamm-Dancoff phonon F^. ISR is generated by the identity: 



Mn - AC, 



(3.33) 



where Mr (t = p, n) denotes the r-particle number operator. Averaging this equation on 
the renormalized ground state and then inserting between the two operators the unity 
operator, one obtains: 



Si = E(0|/^-Mllfc)(l^l/3,|0) 



E ^Oplklbn)(o|AWIU(uK(A')|o)f (vll^lbn') 

p,n;p'n';^ 



it r^'„'Mn\V 



E j{jp\W\ \jn) -J- Up' \W\ \jn')Xk {pn)Xk {pn) .3 

p,n;p'n';k 

Y,{2jp + l){{jp\\a\\jn)fD,{pn). 

p,n 



(3.34) 



Here the index k is an ordering label for the roots of the dispersion equation (3.18). Note 
that Di{pn) is a sum of terms depending exclusively on p and n respectively, and moreover 

ElOpll^lbn)!' = 3, 



ElOnll^lbp)!' = 3. 



(3.35) 
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Taking now into account the result from Eq. ()3.27|) one obtains: 



Sj = 3{N-Z). 



(3.36) 



We may conclude that within the present formalism the Ikeda sum rule is satisfied. A 
peculiar feature of the present result is that ISR is satisfied although the strength of the 
f3~^ transition is vanishing. Therefore, projecting out the components of good particle to- 
tal number from the phonon operator, one obtains a phonon operator which is formally 
of Tamm-Dancoff type and describes a subset of pn excitations obtainable due to the un- 
projected operator. Moreover, ignoring the terms coupling the equations of motion of the 
particle total number conserving operators with those of operators non- conserving the par- 
ticle total number, the two body pp interaction is ruled out. We have seen that assuming 
a consistent decoupling scheme, two independent excitations appear, one describing the (N- 
1,Z+1) system while the other one the (N+1,Z+1) nucleus. The first excitations are suitable 
for describing the rates of the single P~ process while the other type of states are deuteron 
like states and might be populated in a transfer reaction process. Similar formalism can 
be immediately written down for the states in the neighboring (N+1,Z-1), (N-1,Z-1) nu- 
clei obtained in a /?"'" process and by removing a dipole deuteron from the mother nucleus, 
respectively. 

IV. A POSSIBLE IMPROVEMENT OF THE FRPNQRPA APPROACH 

In the previous section we have assumed that the commutators of the operators A"^ and 
At are negligible small. Let us now investigate whether it is possible to define a linear 
combination 



which commutes with both Al^ and A, irrespective of Dp and D„ magnitudes. One finds out 
that a necessary condition is that Dp{k) = Dn{k) which results in having either X{k) = Y{k) 
or X{k) = —Y{k). A similar result is obtained for a linear combination of A^ and A required 
to commute with A^ and A. One may conclude that, rigorously speaking, it is not possible 
to define non-spurious pnQRPA phonon operators as linear superposition of A^ and A which 
are fully decoupled of the phonon operators built up with A^ and A. In this context one 




(4.1) 
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may state that the hnear combination of and A used as phonon operator, in Ref . jisl] . 
which commutes with the particle total number operator is not justified. 
However, it can be checked that there is a linear combination operator 



(4.2) 



corresponding to a specific set of coefficients 



^ p ' p 



1 



1 fu'.-v' 



p ' 



(4.3) 



which commutes with and Ai-^ 



2A(f/2-t/2) V f/nK 

t oTi^ A. a time. Due to this prop- 

erty, it is worth building up a phonon operator out of the independent operators 
{A\Ai-^{—Y^^)]{TZ\i^^,TZi-^{—Y~^). Note the the new proton-neutron operators can also 
be renormalized, 



where 



1 



1 



(4.4) 



D^{pn) = (c^ - z^)Di{pn) + (a^ - b^)D2{pn) + 2{az - bc)D.p + 2(ac - fe)£'„, (4.5) 



so that the new operators satisfy boson like commutation relations: 



ni^{pn),n\^,{pn' 



^p,p'^n,n'^fj,,fi' ■ 



(4.6) 



After some laborious, otherwise elementary, calculations one finds the equations of motion 
for the renormalized operators ^i_^(— )^^^, 7^]^^, 7?.i^_^(— )^^'^: 



H, AUpini] 



H,Ai^^^{pini){-) 



[Tu{pini; pn)A\^{pn) + TuipiUi, pn)Ai-^,{pn){-y ^+ 

pn 

TisipiUi, pn)'JZ\f,{pn) + Ti^{pini]pn)'Ri^^^{pn){-Y'^' , 
[T2i{Pini;pn)A\^{pn) + T22(pini; pri)^i _^(pn)(-)^~''+ 



T2z{pini]pn)Tl\^{pn) + T24,{pini; pn)'JZi-^{pn){- 



H,TZlp{pini)] = J2 [T3i{pini;pn)A\^{pn) + T32{pini;pn)Ai-^,{pn){~Y' 

pn 



H, 7^l^(plr^l)(-) 



TasiPini, pn)n\^{pn) + ^34(^1711; pr^)7^l _^(pn)(- 
Y [T4i{pini; pn)A\^{pn) + T42{pini; pn)Ai,^f,{pn){~y'^+ 



pn 



TAz{pini]pn)TZ\^{pn) + T44^{pini] pn)ni-^{pn){- 



(4.7) 
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Analytical expressions for the coefficients Tik involved in Eq. ()4.7p are given in Appendix B. 
At this stage a phonon operator can be defined: 



rt = ^ [x{k)A\jk) + Z{k)Rl{k) - Y{k)A,,^,m-f-n - W{k)R,,_,{k){-y- 

k=(p,n) 



by requiring that the following equations are obeyed: 



t 



(4.8) 



(4.9) 



The above operator equations provide a system of linear and homogeneous equations for the 
phonon amplitudes: 



A B 

-B* -A* . 









z 




z 




= UJ 




Y 




Y 


\w) 




\w) 



(4.10) 



where the sub-matrices A and B have the following expressions in terms of the matrices T 
involved in the equations of motion written above: 

'Tn f3i\ 

V Ti3 T33 / 



A 



(4.11) 



B 



-^12 -'32 
--^14 -^34- 



(4.12) 



Here Tik denotes the transposed matrix of Tik. 

fik{pini,pn) = Tik{pn,pini). (4.13) 
The second equation ()4.9p yields the normalization equation for the phonon amplitudes: 



Y.[\xm'+\zm'-\Ym'-\wm 



1. 



(4.14) 



The results for the normalization factor D3[pn) are as follows. First one determines the 
boson expansions: 



N, = Y.Gk{p)T,,{k)Tl{k), 

k,fi 

Nn = Y.Gkin)T,^ik)Tlik), 



(4.15) 
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with 

Gk{T) = -(0| [[N^Jr,{k)\ ,rI^(A;)] |0), r = p, n. 
The final result for the average values are: 



(4.16) 



n,k - 



j2 
Jn 



;2 



J 



-E 

p,k 
Jp 



Zlipn) + Wlipn) + Xlipn) + Y^{pn) 
Xlipn) + Y^ipn)) + (Zlipn) + Wl{pn) - Xl{pn) - Y^{pn) 



{Zlipn) + Wl{pn) + Xl{pn) + Y^{pn) 



(4.17) 



X ,Y , Z,W are the amplitudes of the phonon operator when it is expressed in terms of 
the primary operators , A, , B. Their analytical expressions are given in Appendix C. 
From Eq. ()4.17|) one determines the averages (Np) and (Nn) and then by means of (3.7) and 
(4.5) the expressions of Di{pn) and D^lpn) are readily obtained. 

Following a similar procedure as in the previous section, the Ikeda sum rule can be easily 
calculated. The result is: 



Si = J2 hip' {jp\ WWjn) Up' IWl \jn')D^' {pn)D^' (p'n) {Xk{pn)Xk{pn) - Yk{pn)Yk{pn')) , 

pn,p'n';k 

(4.18) 

where the low index k distinguishes the FRpnQRPA solutions. 

If the components Z and W of the FRpnQRPA mode are ignored then the summation 
over k in the above equation would give Spp'Snn' and 



Sj = E(2jp + 1) i{jpMjn)f D.ipn). 



(4.19) 



From this stage on one may follow the path sketched in the previous Section to prove that 



Sj = 3{N-Z). 



(4.20) 



However such a picture is not fully consistent. Indeed, ignoring the amplitudes Z and W 
means to neglect the T13, T14, T23, T24 terms from the equations of motion. To be consistent, 
one has to neglect also the Xi terms form Tu and T12 which would lead to a phonon operator 
to which the two body pp interaction does not contribute at all. Moreover, recall the fact 



19 



that the commutator [A, A] was neglected arguing that factors hke UpVp and UnVn are smalL 
Based on similar arguments, the UV terms from Tn and T12 have also to be left out. Then 
one arrives at the Tamm-Dancoff structure for the phonon operator as we obtained in the 
previous Section. 

Coming back to Ikeda sum rule, it is obvious that the summation over k of 
{Xk{pn)Xk{p'n') — Yk{pn)Yk{p'n')) is an under-unity but close to unity quantity for p = p' 
and n = n', and close to zero when {p',n') 7^ (pn). It results that ISR is slightly underesti- 
mated by the fully renormalized and consistent pnQRPA approach. 



V. SUMMARY AND CONCLUSIONS 



In the previous sections we analyzed two distinct ways of rearranging the two quasipar- 
ticle and quasiparticle density dipole operators inside the FRpnQRPA phonon operators, 
determined by the decoupling hypothesis adopted. In the first scheme we supposed that the 
operators A'^ and A commuting with the particle total number operator, are fully decoupled 
from those which miss this property, i.e. A^,A. Also, we neglected the contribution to 
the equations of motions for the above quoted operators, generated by the pairing inter- 
action which are comparable small with the small terms mentioned before. The result is 
that one obtains two decoupled sets of equations of motion, both of Tamm-Dancoff type. 
One describes a particle-hole dipole excitation in the (N-1,Z+1) nucleus while the other one 
is associated to a deuteron dipole excitation in the neighboring (N+1,Z+1) nucleus. For 
;he first case the ISR is fully satisfied provided the BCS equations are renormalized as well 



IS 



20|. It is proved that, rigorously speaking, to define a phonon operator mixing 



both A'^ and A and keeping at a time the decoupling picture of the operators breaking the 
gauge symmetry associated with the particle total number, is not possible. However it is 
possible to define a linear combination of A'^ , A, A"!" and A denoted by TZ'^ which commutes 
simultaneously with A^ and A. Finally, the phonon operator is built up based on the inde- 
pendent operators A^ , A, TV and 71. We should mention that the original FRpnQRPA is 
lacking this feature. Due to the fact that in the previous publication the effects generated by 
non-commutativity of the component operators of the FRpnQRPA phonon operators were 
thrown out, the results of the present paper constitute a step forward for the renormalization 
formalisms. It is concluded that the ISR is underestimated by the general scheme, but has 
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the advantage that keeps the standard BCS unmodified. One hopes, however, that adding 
a boson expansion (BE) on the top of the FRpnQRPA formahsm, the ISR will be further 
improved. This expectation is supported by the results of a previous work 2l| where a BE 
was performed in terms of the partially renormalized phonon operators (where the scatter- 
ing terms are ignored). The result is determined by the fact that while the renormalized 
pnQRPA underestimates the ISR, the BE has an opposite effect. 



VI. APPENDIX A 



Here we list the exact expressions for the commutation relations of the operators 



{pn) , Ai^ (pn) , AI^ (pn) , Ai^ (pn) : 



A,,{pn),Al,{p'n')\ = J2'i-fClL^,,{-yW{lj,fjnr,jnl) 

f 



A{jpn),A{_Ap'n)i 



6^ 



p,p' 1 



Aif,{pn),Al^,{pn^ 



A,,{pn),Al,{p'n')] = Y^ifC^ 



rrif fi' 



+ w(ijn/v;jpi) 



1 / 1 

fj, mf fi' 



W{ljnfjp';jpl) 



W{ljpfjn';jnl) 



fnif 



^p,p' 



fmf 



Sp,p' 



fmf 



(6.1) 



VII. APPENDIX B 



Here we give the explicit expressions for the coefficients Tik involved in Eq ]4.7l 

Tii{pini;pn) -- 



^Pi ^pi "pi ^^Pi "Pi ^2 t2 



ai-bj 



aj-bi 



^p,pi^n,ni ^ ^"pir!,i;pn 



(1) {biZi — aiCi){bz — ac) + {piCi — aiZi){bc — az) 

X /Ci^pini;pn 



{a^ -b^){al-bl) 



TuipiUi] pn) 



„ biZi — aiCi „ fciCi — Oi^i 

9771 Tj \r ^ ^ ^ ^ A_OTP TJ ^/ 

Pi ^Pi ^Pi 2 l2 ^^ni^ ni ^ n\ 



ai-bi 



al-bl 



^p,pi ^n,ni 



(1) ~ aiCi){bc — az) + (6iCi — aiZi){bz — ac) 



Pini;pn 



{a^ -b^){aj-bl) 

.^3) {biZi - aiCi)a - (6iCi - aiZi)b 
Tis{pini;pn) = -2xicrlJ^,^ 



■,pn 



{a^ -b^){al-bl^ 
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Tu{pini;pn) 
T2i{pini,pn) 
T22{pini,pn) 
T23{pini,pn) = - 
T24{Pini,pn) = - 

T3i{pini,pn) 



(13) -{hzi - aici)b + (&ici - aizi)a 



Pirii :pn 

-Ti2{pini,pn), 
-Tn(Pini,pn), 

-Tii{pini,pn), 
-T*^{pini,pn), 



\Ds{p,n^)\-'/'Dl/\pn) 



a\-hl 



E, 



pi 



v^\)c, + 2u^,vM{cii~hi 



{bz — ac)ai — {be — az)bi 



^Pi,P^ni,n 



- 2xi(T^^^\pini;pn) 



b'M-bl 



T32{pini,pn) 



\D3ip,n,)\-'/'Dl/\pn) 
al-bl 



E, 



{{Ul - Vl)z, + 2U,,V,,a,) {a\ - bl) 
+ {bizi - aiCi) ({Ul - Vl)bi - 2Up,Vp,c,) + {biCi - aiZi) [{U^ - V^ai - 2Up,Vp,Z: 
+ E^, [({V^ - U'Jz, + 2C/„,K,&i) {aj - bj) + {b,z, - a,c,) {{Ul - Vl)b^ + 2[/„,K,^i) 

+ - a,z,) {{Ul - Vl)a, + 2C/„,KiCi); 

(^-,^, , (be — az)ai — (bz — ac)bi 

- 2xia^^^\pini;pny ^ ' ^ ^ ' 

Tsaipifii^pn) 



+ 



^pi,P^ni,n 



(a2-62)(af-6?) 



E.^ 



pi 



{K - KM b'l) - 2UpM(^,z, - b,c,)) 



+ {{Ul - Vl){ai - bl) + 2UnMaic, - b,z^)) 

(3) / ^ acLi + bbi 



- 2xi(J^ '{pini;pn)- 



a^-lP){al-bl 



,2^ ' 



T34{pini,pn) 



+ 2xicr^^\pini;pn) 



bai + abi 



[2Ep^Up^Vp^{biZi - aiCi) + 2£;„^C/„i (oi^i - biCi)] 



{a^-b-'){al-biy 
T4iipini,pn) = -r32(pini,pn), 

T42{pini,pn) = -T^^{pini,pn), 

Ti3{pini,pn) = -T^i{pini,pn), 

T44{pini,pn) = -T^2(j)ini,pn). 



(7.2) 



(7.1) 
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The expressions for a, b, c, z have been defined before by Eq.(4.3). The coefficients ai,bi,ci, z\ 
are obtained from Eq.(4.3) by replacing the indices p and n by pi and ni, respectively. The 
index * stands for the complex conjugation. Several notations have been used: 

cr^^^(pini,pn) = ep^„Ji:)3(pini)|^/Vp^„^,p„|i:'3(pn)|^/^ 

cr(^^)(pini,p?2) = I5i/^(pini)ap^„j,p„|D3(pn)|^/^ 

= ep^„jL'3(pini)|^/Vpi„i,p„Di^^(pn), 
_ I^3(pini) 

"^"^ ii^3(pim)r ^'-^^ 

VIII. APPENDIX C 

Here we give the explicit expressions for the bar amplitudes involved in the equations for 
the quasiparticle number operators averages ()4.17|) . The results are: 

1 1 

^fc(pn) = , {UpVnXkipn) - UnVpYk{pnj) 



Di{pn) ^J\D3{pn)\ 
[Zkipn) {allpUn - bVpVn + cUpVn + zUnVp) 

- Wk{pn) {-aVpVn + bUpUn + cUnVp + zUpVn)] , 

-%{pn) = , ^ {UnVpXk{pn)~UpVMpn)) ^ 



Di{pn) ^\D;>,{pn)\ 
[Zkipn) {-aVpVn + bUpUn + cUnVp + zUpVn) 

- Wk{pn){aUpUn-bVpVn + cUpVn + zUnVp)], 

Zkipn) = (UpUnXkipn) + VnVpYk{pn)) + 



Di{pn) ^J\D-i{pn)\ 
[Zkipn) {-aUpVn - bVpUn + cUpUn - zVnVp) 

- Wk{pn) {-aVpUn - bUpVn - cVnVp + zUpUn)] , 

-Wk{pn) = , ^ {-VnVpXk{pn)-UpUnYk{pn)) 



Di{pn) ^[D3{pn)[ 
[Zkipn) {-aVpUn - bUpVn - cVnVp + zUpUn) 

Wk{pn) {-aUpVn - bVpUn + cUpUn - zVnVp)] . (8.1) 
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